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X ■ Abstract 

Baxter's Q-operator is generally believed to be the most powerful tool for the exact diago- 
nalization of integrable models. Curiously, it has hitherto not yet been properly constructed 
in the simplest such system, the compact spin-- Heisenberg-Bethe XXX spin chain. Here 
we attempt to fill this gap and show how two linearly independent operatorial solutions to 
Baxter's TQ equation may be constructed as commuting transfer matrices if a twist field 
is present. The latter are obtained by tracing over infinitely many oscillator states living in 
the auxiliary channel of an associated monodromy matrix. We furthermore compare and 
differentiate our approach to earlier articles addressing the problem of the construction of 
the Q-operator for the XXX chain. Finally we speculate on the importance of Q-operators 
for the physical interpretation of recent proposals for the Y-system of AdS/CFT. 



1 Motivation 

Recently there was much progress with integrability in planar four-dimensional gauge the- 
ories [U EJ [3] and AdS/CFT jH |5]. At variance with the long-held belief that quantum 
integrability is confined to low-dimensional systems, the asymptotic Bethe Ansatz solu- 
tion of planar M = 4 gauge theory was conjectured by a combination of rigorous results 
and assumptions [5J. Further conjectures and arcane techniques led to first proposals [6] 
for the full, non-asymptotic spectrum of this gauge theory, as well as its dual superstring 
theory on the space AdSs x S 5 . The proper interpretation and even the veracity of these 
latter proposals remains hotly debated, and in any case there is currently no trace of the 
theoretical underpinnings of the employed experimental mathematics. The final equations 
obtained take the form of an infinite system of integro-difference equations, based on a 
conjectured thermodynamic Bethe Ansatz, which involve an infinite set of "T-functions" . 
These may in turn be rewritten for a sometimes more convenient set of "Y- functions" . 
In the large volume limit the "T-functions" are to be interpreted as transfer matrices of 
an infinite number of excitations, most of them corresponding to bound states. However, 
some of the T-functions involve excitations believed to be elementary (mirror-magnons) , 
while others turn into what looks like the eigenvalues of some Baxter Q-operators. 

It is currently totally unclear which operators T and Q possess all these eigenvalues 
T and Q. In other words, even if the proposals |6J are made more precise, and turn 
out to be correct, the question "What has been diagonalized?" will remain unanswered. 
One clearly would like to construct the associated operators, and prove that their mutual 
commutativity is based on an underlying Yang-Baxter symmetry. 

One very curious feature of AdS/CFT integrability is that the underlying integrable 
system is both a quantum sigma-model "living" on a smooth continuous two-dimensional 
worldsheet, and at the same time a certain long-range spin chain defined on a discrete 
lattice. The two pictures are related by a continuous coupling constant, and one cannot 
say that the continuum sigma model is obtained from a discrete spin chain by a continuum 
limit. The sigma model is also a spin chain. 

For short-range spin chains, the construction of transfer matrix operators T is rather 
well understood, see below. However, the general principles of construction of the Baxter 
Q-operators, originally introduced in his seminal paper on the 8- vertex model [7], are less 
well understood, and appear to be less systematic, even though the subject was intensively 
studied for the past twenty years. An important step toward a general algebraic theory of 
the Q-operators, particularly relevant for this paper, was made in [SI [9] devoted to confor- 
mal field theory where these operators were constructed as traces of certain monodromy 
matrices, associated with infinite dimensional representations of the g-deformed harmonic 
oscillator algebra. More generally this method applies to any model with U q (sl(2)) symme- 
try and it was further generalized for the case of some higher-rank quantized algebras and 
super-algebras [TQl [Hi [12] . However despite all these considerations of rather complicated 
models with "g-deformed" symmetries, it appears that there is still no completely explicit 
construction in the literature for the Q-operator of the compact XXX chain - - the first 
spin chain ever solved by Bethe Ansatz [T3]. We therefore decided to fill this gap, as a 



very first necessary step to vigorously address the much more involved case of AdS/CFT 
integrability. Excitingly, we will find that we need infinite dimensional representations to 
carry out this construction. These are therefore needed to properly understand the in- 
tegrable structure of spin chains with finite dimensional quantum space, a feature which 
then puts spin chains on a par with integrable sigma models. We believe that this makes 
the appearance of spin chains in the strongly coupled limit of the AdS/CFT sigma model 
much less surprising. 

2 Brief Review of the Spin-^ XXX Heisenberg Chain 

The one- dimensional Heisenberg spin chain Hamiltonian reads 

L 



H = 4^Vi-Si.£, +1 ] with S=^a, 

1=1 ^ ' 



(2.1) 



where a are the three Pauli matrices, i.e. S is the spin-| representation of su(2). This 
Hamiltonian acts on the L-fold tensor product 

V = C 2 ®C 2 ®---®C 2 , (2.2) 

L— times 

which will be called quantum space throughout the paper. The spin operator Si acts on 
the Z-th component of the quantum space, and it is clear from (12. ip that we need to specify 
the meaning of Sl+i- Periodic boundary conditions are imposed by defining 



>L+1 



Si . (2.3) 



-L 



This Hamiltonian may be rewritten as H = 2 Xw=i (Ji,i+i ~ ^M+i) where \i+i and Pz,z+i 
are the identity and the spin permutation operators on adjacent sites (1,1 + 1) of the chain of 
length L, respectively. It appears in M = 4 Yang-Mills theory in the scalar field subsector, 
where S G su(2) C su(4) C psu(2, 2|4), as the one-loop approximation of the conformal 
dilatation generator D e su(2, 2) C psu(2, 2|4) 

B = L + g 2 U + 0(g 4 ). (2.4) 

Here g 2 is related to the 't Hooft coupling constant A by g 2 = j^z- Note that the Hamil- 
tonian (12.11) with boundary conditions (12. 3 p is rotationally invariant, i.e. [H,«S] = 0. 

It is well known that Hans Bethe discovered in 1931 a system of algebraic equations 
which yield the exact spectrum of H. This was obtained after making an Ansatz for 
the wave function now carrying his name JT3]. This so-called coordinate Bethe Ansatz 
interprets the state t • • • t with energy eigenvalue E = as the vacuum. Each up-spin | is 
an unoccupied lattice site, and each down-spin J, is interpreted as a lattice particle, termed 



magnon, carrying lattice momentum p. After introducing a rapidity u for each magnon, 
where 

e ip = ^±i ^ u =l cot l (2.5) 

u-\ 2 2' K J 

Bethe's solution for the eigenvalues E of H in the (conserved) sector of M magnons reads 

- 1 

^ = 2 E^-T' ( 2 - 6 ) 

fc=l Uk + 4 

where the M Bethe roots u^ have to satisfy the Bethe equations 

, i r M 

Uk+o\ TT u k~ Uj + 1 



n fc ~ ^ • (2.7) 



Wfe - | / ^ U k - Uj - I 

Uk 

The eigenvalue U of the lattice translation operator U, which shifts a given spin configu- 
ration by one lattice site, is given by 

M _ {_ 

There has been a longstanding controversy, starting from a discussion in Bethe's original 
paper, as to whether the obtained spectrum is complete, which still attracts much attention; 
here we give only a small sample of references [T^l [T5l [T6l [T7j . There are (at least) three 
distinct subtleties. The first is, that the map f)2.5p leads to infinite rapidities u = oo for 
zero-momentum magnons where p = 0. Of course this is consistent with both the expression 
for the energy (12.61) and the Bethe equations (12. 7p . but complicates the proper counting 
of states. The effect may be traced to the su(2) invariance of H: The highest weight 
state of each energy multiplet corresponds to a solution with only finite rapidities. The 
descendents of this state are obtained by applying the global su(2) lowering operator. Each 
application places a further p = magnon into the chain, corresponding to a completely 
symmetrized insertion. A second subtlety is that it is a priori not allowed to increase M 
beyond half-filling, i.e. such that M > L/2 (all states must then be descendents). Spurious 
states with finite rapidities are nevertheless found. In fact, one "experimentally" finds that 
solving the equations 

C*±X\ = TT Uk -^ + i (2.9) 

\u k -^J £1 Uk - u .- % 

Uk 

for L — M + 1 roots u^ and plugging the solution into 

L-M+l 

E = 2 J2 tfTXT ( 2 -!0) 

fe=l U k + 4 



yields energies identical to the ones of highest weight states with magnon number M. 
However, the root distribution Uk = Uk{ot) depends on one arbitrary complex parameter a, 
while the energy does not. These "beyond the equation solutions" were discussed in detail 
in |18j . The third subtlety, already noted in [13], is that some momenta appear in pairs 
of the form p = po ± i oo, which leads via (12. 5p to Bethe roots at u = ±|. The physical 
picture here is that two magnons form an infinitely tight bound statqj, they are "stuck 
together" . However, they of course contribute only a finite amount of energy, so (I2.6P has 
to be interpreted with great care. 

It is also quite well known that all these subtleties may be resolved by replacing H by 
a "twisted" Hamiltonian H^, where can be interpreted as an (imaginary) "horizontal 
field" in condensed matter parlance, or alternatively as a magnetic flux passing through 
the chain looped into a circle 



H * = 4 E i\ ~ S ? 5 w " \ e%i S ? ^+1 - ~ e-** ST 6J.J , (2.11) 

i=i 

with iS z = Si ±iSf. The Bethe Ansatz still works with minor modifications. The formula 
for the energy (12.61) and the relations ( 12. 5 p remain unaffected, but the Bethe equations 
( 12. 7p are modified to 

g*±lV> = FT* ~ "' • (2-12) 

u k - | / ** Uk -Uj-i 

This tiny modification resolves all the difficulties we just discussed (for generic values of 
the twist). Firstly one now finds all 2 L states of the length L spin chain, such that all 
rapidities u k are finite. Secondly, there is no "beyond the equator problem" anymore and 
M can range over M e {0, 1, . . . ,L — 1,L}. In fact one can now rigorously derive the 
second Bethe Ansatz, which treats the up-spins as "dual" magnons in the vacuum of the 
down-spins. It reads 

- . i T L-M _ 

' UJ1±A ) V* = TT Uk " Uj — , (2-i3) 

Uk - | ' J7 u k -Uj — i 

Uk 
L-M 



E = 2 E^TT> (2-14) 



k=1 K + 5 



and in contradistinction to (12. 91) . (12. 101) there are only L — M dual magnons, as expected. 
In fact, the set of equations (12.121) and (12.131) are completely equivalent, each set yields 
the full spectrum of all 2 L states. Finally there are no more degenerate roots at u — ±| 
with a priori indefinite energy ( I2.6P or ( I2.14p . The price to pay is that su(2) invariance of 
the spectrum is broken: All multiplets split ud_|. However, we can think of as a small 



1 This happens first for one of the two singlet states of the L = 4 spin chain. 

2 A spin s multiplet where the magnetic quantum numbers are m = —(2s + 1), . . . , (2s + 1) splits up 
such that there still is a degeneracy between any two states whose m differs by a sign flip. 



regulator which may always be removed where physically sensible. Note that such twists 
are also natural from the point of view of the AdS/CFT correspondence. They appear in 
the scalar sector of the integrable one-loop dilatation operator of the /3-deformed twisted 
Af = 4 gauge theory [19] . 

The magnetic flux <fi may be distributed in many possible ways. For instance, as 
concerns the energy spectrum of the spin chain, it is equivalent to use instead of (12. lip a 
linearly transformed Hamiltonian, 

H^ = C((j>/L) H^ C((j)/L)~\ C(a) = e iLaS 'i <g> e *(^-i)««!-i <g, . . . <g, e iaS " , (2.15) 

which is given by the original formula (12. ip . but with "twisted" boundary conditions 

Sl +1 :=Sl St +1 :=e^Sf. (2.16) 

There is an interesting way to reformulate the Bethe equations ( 12.12)) and f !2. 13|> with 
the help of Baxter polynomials defined for each eigenstate by 

M L-AI 

A_(tt):=IJ(«-u fc ), A+(u):= Y[(u-u k ). (2.17) 

fe=i fc=i 

Then 

( f-j e ±l * = — 7^7 if u = u k or u = u kl respectively. (2.18) 

Vm-|/ A T (u — i) 

On infers that the following equations must hold for each eigenstate and for all mgC: 

T(u) A T (u) = e ±l * (u + l -j A^u - i) + e** (u - % - J A T (u + i). (2.19) 

The reason is that the r.h.s. is a polynomial with roots at u = u k or u = u k , respectively, so 
it must be proportional to A^(u). We also see that the latter function must be multiplied 
by another polynomial T±(u) of degree L. One experimentally finds that this further 
polynomial does not depend on ±, i.e. T±(u) = T{u), a fact which will be proven in the 
next chapter. So we have from (I2.19P and (12.171) 

6 L 

T(u) = 2cos| Y[(u-w k ), (2.20) 

fc=i 

where the w k are the remaining L roots of the r.h.s. of (I2.19p . 

Equation (I2.19P can be regarded as a second order difference equation for an unknown 
function A(u), which has two linearly independent solutions A±(u). To see this more 
clearly, it is convenient to define the Baxter functions 

Q ± (u) := e ±u % A ± {u) , (2.21) 



such that Q± (u) are indeed two linearly independent solutions^ of the difference equation 

T{u)Q{u)= (« + £) Q(u-i)+(u-^J Q{u + i) with Q(u) = Q±{u). (2.22) 

This is Baxter's famous TQ equation for the twisted Heisenberg magnet. As we can see, the 
twist has actually disappeared from the equation, and is entirely encoded in the analytic 
Ansatz fl2.2ip . fl2.17p . fl2.20p for the solution. Note that the Baxter functions at nonzero 
twist <p are not polynomials. 

Baxter derived this equation, which holds for all eigenvalues of the commuting T- 
and Q-matrices on the operatorial level, in his original solution on the "zero-field" 8- 
vertex model [TJ, which also contains the solution of the "zero- field" XYZ spin chain |2U] . 
Although it seems to be possible to take a limit of his results and apply them to the 
untwisted XXX model, this would then only apply to the "zero-field" case = 0. 

One purpose of this article is to consider the XXX chain with an arbitrary non-zero 
field 7^ and provide an independent construction of the operators T(u) and Q±(w) 
satisfying f)2.22p such that all their eigenvalues are of the form spelled out in f!2.2ip , ( I2.17P 
and (I2.20p . Our construction of the Q-operators is conceptually very similar to that of 
[SJ E] in c < 1 conformal field theory, based on the Yang-Baxter equation with U q (sl(2)) 
symmetry. However, here we provide a self-contained and separate consideration of the 
XXX model, rather than attempting to take the q 2 — > 1 limit of the relevant results of 

Eli- 

In the case of the operator T(u) Baxter's construction actually immediately applies. 
His methodology was subsequently developed and systematized by the inverse scattering 
methodology of the Leningrad school of mathematical physics. Here we will just briefly 
summarize the construction, referring for all further details to the authoritative presenta- 
tion [2T] . One constructs this operator, termed transfer matrix T(w), as the trace over a 
certain monodromy matrix M.(u). The latter is in turn built from a "generating object", 
the local quantum Lax operator 

u + iSf i Si 



£ <M=I iSf' u-iSfh (2 ' 23 > 



which is a 2 x 2 matrix in some auxiliary space C 2 with site-/ spin operator (as in (12. ip ' 
valued matrix elements. The monodromy matrix is then built aaj 

M(u)=[ J ' °.i | ■ C L (u) ■ Cl-i(u) ■ . . . ■ C 3 (u) ■ Ci(u) . (2 2 1! 

e 



3 The most general formal solution of (|2.22[) is then a linear superposition of the form Q(u) = 
c + (u) Q+{u) + c_(n) Q-(u). However, unlike the theory of second order differential (as opposed to differ- 
ence) equations is that the "constants" c±(u) could a priori be any functions of u with period i. 

4 Hcre • denotes 2x2 matrix multiplication in the two-dimensional auxiliary space. The entries of this 
2x2 matrix act on (|2.2|) . 



In view of (12. 2p it is an operator acting on the tensor product of the auxiliary space and 
the quantum space. Finally, by taking the trace over the two-dimensional auxiliary space, 

T(u) = TiM(u), (2.25) 

one constructs the transfer matrix as an operator on the quantum space (12. 2p . It is easy 
to show, see [2T], that at the special point u — \ the transfer matrix becomes proportional 
to the lattice shift operator (cf. discussion around (I2.8P ). multiplied by a diagonal matrix, 
i.e. one has 

U = r L T(f)e- i05 " . (2.26) 

The Hamiltonian ( I2.15P is then obtained from the expansion of the transfer matrix in the 
vicinity of the point u = | , 



H^ = 2L-2z — logT(w) 

du 



(2.27) 



Finally, because of the underlying Yang-Baxter symmetry, to be discussed below, one can 
show that transfer matrices with different values of spectral parameters form a commuting 
family 

[T(«),T(«')] = 0, (2.28) 

which also contains the Hamiltonian (I2.15p . obtained from (12.11) by twisting the boundary 
conditions (I2.16p . However, this fact by itself does not directly lead to the solution of 
the model. One procedure is to apply the algebraic Bethe Ansatz as explained in [2~T] . 
Another, to be developed below, is to also construct the operators Q±(w) as traces over 
some monodromy matrices, and to derive the operator version of (I2.22p 

T(u) Q±(u) = (u+ \) L Q ± (u - i) + (u - f) L Q ± (u + i) , (2.29) 

along with a proof that the analyticity of the eigenvalues of the therein appearing operators 
is given by ( l2.2ip .( l2~TT7]UI2.20p . Then no Bethe Ansatz is required, and the twisted Bethe 
equations (I2.12p or ( I2.13P immediately follow (for more details on this logic, please see 
section 13.51 below) . 

A very interesting issue is the (f> — > limit of the operators Q±(w) appearing in (I2.29p . 
In this limit, the broken su(2) invariance of the spin chain is recovered. Then the majority 
of the eigenvalues of Q±(w) turn into descendents of su(2) highest weight states. From 
our discussion following the untwisted Bethe equations (12. 7p these will have Bethe roots at 
u = oo, so the operators Q±(u) are expected to diverge, in stark contrast to the eigenvalues 
of T(m), which are perfectly finite and smooth in the zero-twist limit. On the other hand, 
we can of course retrace the steps leading to the existence of the operator equations ( I2.29p 
by making, in analogy with ( I2.17p . a polynomial Ansatz for linearly independent Baxter 
functions^ 

M L-M+l 

Q(u)~H(u-u k ), P(u)~ J] (u-u k ). (2.30) 

k=l k=l 



5 Note that the exponential factors in (12.211) disappear at <p = 0. We have used the notation ~ since 
the proper normalization of these functions will turn out to be quite subtle. 



This empirical way had already been explored in [T3], but it was erroneously (as we shall 
see) stated that P(u) and Q(u) are completely unrelated to the eigenvalues of Q±. By 
using the untwisted Bethe ( 12 .7)1 and dual Bethe f 12 . 9 j) equations, the latter being the reason 
for the power L — M + 1 in the polynomial P(u), we conclude that there must be finite 
operators Q(u), P(w) which should satisfy the same Baxter equation (12.291) as Q±(w), such 
that however (1) their eigenvalue spectrum is su(2) invariant, and (2) their eigenvalues are 
indeed of the form in (I2.30p . It is clear that finding Q(m),P(m) from Q±(u) must indeed 
be quite non-trivial, and must involve some kind or "renormalization" of these divergent 
operators. It is furthermore a priori quite mysterious how the extra (L — M + l)-th root 
appears in P(u) in (12.301) as compared to the L — M roots of A+(u) in ( I2.17p . These puzzles 
will be resolved in chapter HI In particular, we shall find the resolution to be intimately 
connected to the exponential factors in (12.211) . 

Our article is organized as follows. In the ensuing chapter [3] we will construct the 
operators Q± as the trace over an appropriate monodromy matrix. This will require the 
introduction of two copies of infinite oscillator Fock spaces, despite the fact that we are 
dealing with a finite dimensional spin chain carrying finite dimensional representations of 
su(2). In the next chapter H] we will study the very subtle <p — V limit of our construction, 
leading to a one-parameter family of linearly independent operators Q(w), P{u) by "renor- 
malizing" the previously obtained operators Q±. We will also perform some numerical 
study of their spectrum for a few cases, illustrating an interesting pattern between their 
respective eigenvalue root distributions. In the following chapter |5] we briefly discuss how 
our construction and result differ from a large earlier, complementary literature on the sub- 
ject. We end in E] with a brief list of the many open problems and potential applications 
related to our result. 

Before proceeding with the construction of the operators Q±, let us however change 
notation by 'Wick-rotating" the spectral parameter u to 

z:=-iu. (2.31) 

It is true that the ^-convention of this review chapter [2] is the most widely accepted notation 
in much of the Bethe Ansatz literature, and nearly all of the literature on AdS/CFT 
integrability. But the z-convention (I2.3ip used for the rest of this paper (with the exception 
of the numerical work presented in section 14.51 and appendix |Fj as well as the examples of 
small chain lengths in appendix |B]) prevents all further derivations and manipulations of 
functional equations from being cluttered by factors of i. This will then turn ( I2.29P intd_| 

T(z) Q±(*) = (z + \f Q ± (z - 1) + (z - \) L Q±(z + 1) . (2.32) 



6 For simplicity, and with slight abuse of notation, we will not use new symbols for the various operators. 
E.g. the transfer matrix in (|2.29p is i L times the transfer matrix in (|2.32[) . It should be straightforward 
to return to the notation with the spectral parameter u whenever needed by using (|2.31[) . 



3 Construction of Q± as Transfer Matrices 

3.1 Yang-Baxter equation and commuting transfer matrices 

As explained above the Hamiltonian (12. ip . with twisted boundary conditions (I2.16p . is 
generated by the spectral parameter-dependent transfer matrices (I2.25p . which form a 
commuting family, cf. (I2.28p . Here we want to construct further important transfer matrices 
which nevertheless belong to the same family. To this end one needs to study the possible 
solutions of the Yang-Baxter equation 

R(x - y) (h v {x) <8> 1) (1 ® hv(y)) = (l ® W(y)) (h v (x) ® l) R{x - y), (3.1) 

where R(z) is the rational 4x4 i?-matrix, 

R(z): C 2 ®C 2 ^C 2 ®C 2 , R{z) = z + R, (3.2) 

and the L-operator Ly (z) is a 2 x 2 matrix, acting in the quantum space of a single spin-|, 

Its matrix elements are operator- valued functions of the variable z, acting in an auxiliary 
vector space V. The .R-matrix acts in a direct product of two-dimensional spaces C 2 ® C 2 
and the operator P permutes the factors in this product. Note that the .R-matrix (13. 2p is 
GL(2)-invariant 

R(z) = (G®G)R(z)(G®G)-\ GgGL(2), (3.4) 

where G is any non-degenerate 2x2 matrix. 

The solutions of (13. ip which we will use in this paper are rather simple. The first one 
is a well-known generalization of (I2.23p . 

r3 



z + 3 A J 
z 



H*)= T+ /- T3 = 2l + 2 E^ jfc ' ( 3 - 5 ) 



k=i 



where J 1 * 1 = J 1 ± iJ 2 and J 3 are the generators of the sl(2) algebra 

[J 3 ,J ± ]=±J ± , [J + ,J"]=2J 3 . (3.6) 

In the second form of L(z) in (13. 5p we used the 2x2 spin operators S h appearing in (12.11) . 
Note that for the L-operator (I3.5P equation (13.11) holds on the algebraic level in virtue of the 
commutation relations (13. 6p . To obtain a specific solution one needs to choose a particular 
representation of these commutation relations. For further reference define the highest 
weight representions of sl(2), with highest weight vector t> , defined by the conditions 

3 + v = 0, J 3 v =jv , (3.7) 



where j is the spin. The (2j + l)-dimensional representations with integer or half- integer 
spin, i.e. 2j e Z> , will be denoted by ttj, while infinite-dimensional representations with 
arbitrary complex spin 2j 6 C with be denoted as Kj. 

For each solution of (13.11) one can define a transfer matrix 

T v (z) = Tr v (p Ly(z)®L v (z)(g)---®Lv ) , (3.8) 

L— times 

where the tensor product is taken with respect to the quantum spaces C 2 , while the operator 
product and the trace is taken with respect to the auxiliary space V. The boundary twist 
operator D, specially suited for our purposes, is defined by 

[D, Ln(z)] = [B,L 22 (z)} = 0, DLi 2 (z) = e^Li 2 (z)D, PL 21 ( z ) = e^L 21 (z)D, 

(3-9) 
where Ly(z), i,j = 1, 2, denotes matrix elements of Lj/(z). Note that D is an operator in 



1 



the auxiliary space; it acts trivially in quantum spacej. For the solution (13. 5 j) one easily 
finds 

D = e i0j3 . (3.10) 

Substituting the last expression together with (13. 5p into the definition (13. 8p and taking the 
trace over the standard (2j + l)-dimensional representations ttj, one obtains an infinite set 
of transfer matrices 

T J (z) = Tr 7Tj (M(z)), (3.11) 

built from the monodromy matrices 

M(z) = e^ 3 '" L{z) ® L{z) ® ■ ■ ■ ® L{z) , 2j e Z> , (3.12) 



L— times 



and labelled by the value of the spin j = 0, ~, 1, |, . . . , oo. Note in particular thatj 

T(z) = Ti(z) (3.13) 



2 



coincides, up to a trivial rescaling of the argument z = —iu, with the transfer matrix 
T(u) defined in (I2.25p . The transfer matrices (I3.1ip depend on the spectral variable z and 
(implicitly) on the twist parameter <fi. 

Standard arguments based on the Yang-Baxter equation immediately imply that the 
operators Tj(V) belong to a commuting family, since one derives from (13.11) in generalization 
of fl2T28|) thafl 

[T(z),T J (z')] = 0, 2jeZ> . (3.14) 



Actually, all considerations presented in this paper are valid for a more general case, when the operator 
D acts non-trivially in the quantum space. The only additional condition is that it should commute with 
the diagonal elements A4jj(z), j = 1, 2 of the monodromy matrix (|3.12p . For example the twist could be 
of the form <fi — a + /3Sf ot , where a and /3 are c- numbers. This type of twist originally arose in |22) in 
considerations of commuting T-operators in CFT. 

8 As opposed to (|2.23[) and (|2.24[) . the 2x2 matrices L(z) in (|3. 12|) act on the corresponding copy of 
C 2 in the quantum space (J2.2I) . 

9 One can prove that the more general relations [Tj(z),Tj/(z 1 )] = with 2j,2j' £ Z>o also hold. 
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The presence of the boundary twist f l3.10|) in the definition (13. lip does not affect the 
commutativity thanks to the properties (13.91) and (13.41) . 
In similarity to ( 13. lip we may define transfer matrices 

Tf(z)=Tr v +(M(z)), 2jeC, (3.15) 

where Ai(z) is the same as in (13.121) . but now the trace is taken over an infinite- dimensional 
representation 10 ! tc^ with an arbitrary, possibly complex spin j. This representation is 
spanned by the vectors {vk}^L , with the following action of the generators 

J 3 v k = (j - k) v k , J~v k =v k+1 , J + v k = k (2j - k + 1) Vfc-i • (3.16) 

The convergence of the trace in (I3.15P requires one to assume that Im0 < 0. For a generic 
value of j G C the representation ( I3.16P is irreducible. However when j takes non-negative 
(half) integer values 2j e Z> this representation becomes reducible. The matrices 7r+( J + ), 
tt^(J~) and 7r^(J 3 ) then acquire a block-triangular form with two diagonal blocks. One of 
these is finite-dimensional, being equivalent to the (2j + l)-dimensional representation ttj, 
and the other one is infinite-dimensional, and coincides with the highest weight represen- 
tation 7ri-_i- Hence on the level of traces one easily obtains 11 ! 

Tfc) = Tj(z) - Tt^z) , 2j G Z> . (3.17) 

3.2 Functional relations 

Our immediate aim is to study various algebraic properties of the transfer matrices ( 13. lip 
and ( I3.15p . and, in particular, to derive all functional equations they satisfy. In the context 
of integrable models, related to the quantized Kac-Moody algebra U q (sl(2)), this problem 
has been previously solved in [HI [9] . Here we will apply a similar approach. The detailed 
considerations of [H [9] are devoted to conformal field theory and cannot be straightfor- 
wardly used for the lattice XXX-model. Moreover, the latter model is related to a rather 
subtle limit q 2 — > 1 in the relevant g-deformed constructions of [SI IS], and requires ad- 
ditional considerations. Nevertheless one should expect that the general structure of the 
functional relations, independent of the value of the deformation parameter q, must remain 
intact in the q 2 — > 1 limit by continuity arguments. In particular, one should expect from 
Eq.(2.45) of [S\ that the transfer matrix T^(^), defined in (13.151) above, factorizes into a 
product 

f{<f>) T+(z) = Q + (z + j + i) Q_(z -j-\), where f{<f>) = 2i sin t , (3.18) 

of two Baxter Q-operators Q±(z), which should satisfy 

[Tjiz), Q ± (z')} = 0, [Q + (z) , Q_(^)] = , [Q ± (*) , Q±(^')] = . (3.19) 



10 In general this representation is not unitary, but this is not relevant to our present construction. 
u By the methods of appendix [Cl the relation p. 171) may also be analytically continued to complex j, 
such that Tj(z) stays finite in the <f> — > limit. 
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They clearly extend the commuting family ( 13.14=)) . and should satisfy the TQ-equation 
(j2.32p . Below we will prove that this is indeed true in our case and give the corresponding 
definitions of the operators Q±(z). 

Before going into this proof let us demonstrate that the relation (13.181) alone leads to 
a simple derivation of all functional relations, involving various "fusion" transfer matrices 
Tj(z) and Q-operators [SJ [H]- For this reason Eq. (13.181) can be regarded as a universal 
fusion relation — once it is derived, no further algebraic work is required. 

Substituting (13.181) into (13.17)) one obtaino 

f(<j>)T j (z) = Q+(z+j^)Q4z-j-±)-Q4z+j+±)Q+(z-j-l), 2jeZ> . (3.20) 

For j = (trivial representation) one clearly has 

T (z) = z L I, (3.21) 

where I is the identity operator on the quantum space, which hereafter will be omitted. In 
this case (13.201) reduces to 

Q+(^ + §) Q_(2 - |) - Q.(z + \) Q + {z -\) = /(0) z L . (3.22) 

Using the last relation together with the expression (13.20)) for Ti(z) = T(z) one imme- 
diately derives the TQ-equation ( I2.32() . Note that the relation ( 13.22)) can be regarded as 
the (quantum) Wronskian relation for the second order difference equation ( I2.32p . ensuring 
linear independence of the two solutions Q + and Q_. Furthermore, ( I3.20p clearly indicates 
that the operators Q± should be considered to be more fundamental than the transfer 
matrices, since the latter are a quadratic superposition of the former. 

In fact, we can consider (13.20)) as the most fundamental fusion relation of the model, 
from which all other relations follow. The mechanism by which this happens is quite 
simple. Let us define, for any commuting quantity Q^, 

Tab = Q [j4 Qsj j (3.23) 

where the square brackets indicate antisymmetrization and the indices A, B denote very 
generally a collective set of discrete indices and/or continuous variables. The latter are 
related to linear combinations of the the spectral parameter z and some representation 
labels such as the spin j. We immediately infer the identities 

^Qc] = ' T [A bT c] d = 0. (3.24) 

These two types of equations are a compact way to respectively write Baxter's equation 
( 12.32)) as well as all fusion relations ( 13.25)) such as the ones of |24) : 

T J (z + l)T J (z-l) = (z + j + 1 1 ) L (z-j-l) L + T J+l2 (z)T J _ l2 (z), j = 0, |, 1, . . . 

(3.25) 
A generalization of these relations appears in appendix [A] The same kind of generalized 
fusion relations can be found in 



12 This equation appeared for the XXX model in [23]. However, there the Q-operators were not completely 
explicit cit. "the operators Qjj(0) are not easily determined and we are missing at the moment concrete 
expressions for them." 
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3.3 Factorization of the /^-operator 

As noted before, we want to construct the Q-operators as transfer matrices (13. 8p built from 
suitable //-operators, solving the Yang-Baxter equation ( I3.1J1 . We shall soon see that the 
required L-operators indeed exist. They are easily obtained via some special reductions of 
( 13. 5p . In preparation of this calculation we need to recall some well-known realizations of 
the sl(2) commutation relations (I3.6P in terms of the harmonic oscillator algebra 

H: [h,a ± ] = ±a ± , [a~,a + ] = l, h = a+a + \ . (3.26) 

The Fock representation T of this algebra is spanned on the vectors {vk}^L Q , 

7: a + v k = v k+ i, a.~Vk = kvk-i, hv k = (k + \) v k . (3.27) 

The value of the quadratic Casimir operator of the algebra s ((2) 

C 2 = J 2 = (J 3 ) 2 + I (J+ J- + J- J+) (3.28) 

for the highest weight representations (13. 7p is given by 

7r,(C 2 ) = 7r/(C 2 ) = j(j + 1) ■ (3.29) 

Below we will use the fact that for the infinite-dimensional representation 7r+ the sl(2)- 
generators 

J? = tt+(J°), a = {3,+,-}, (3.30) 

can be realized through the oscillator algebra (j3.26p 

jT = a+, J+ = (2j-a+a-)a-, Jf = j - a + a~ , (3.31) 

or alternatively as 

jT = a + (2j-a+a-) , J+ = a" , J 3 = j - a+ a" , (3.32) 

where the operators a + , a - are taken in the Fock representation (13.271) . These realizations 
are commonly known as Holstein-Primakoff representations; they can be readily verified. 
Indeed, having in mind the Fock representation (I3.27p . one can easily see that the formulae 
(I3.3ip are essentially verbatim transcriptions of the matrix elements (I3.16P in the oscillator 
notations. The second realization (J3.32J) is obtained from (13.311) by a simple similarity 
transformation 13 !, and, therefore gives the same trace in (13.15jl 1 'i. 



13 One can also take 

J7=a+(2. j --a+a-) 7 , J+ = (2j - a+a-) 1_7 a- , J| = j - a+ a" . 

with arbitrary < 7 < 1. Here we only use the values 7 = 0, 1. 

14 This is also true for 2j € Z>o, when (|3.31[) and (|3.32|l become reducible and on the level of traces both 
lead to the same formula (I3.17[) 
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Due to the invariance of the R- matrix (13.4)) . the solutions of (13.1 p are defined up to 
linear transformations 

h v (z) -> F h v (z) G, F,Ge GL(2) , (3.33) 

where F and G are arbitrary non-degenerate 2x2 matrices. In other words, the transfor- 
mation (13.331) does not affect the validity of the Yang-Baxter equation (13. ip . 

Using the realization (13.311) . let us write explicitly the L-operator (13. 5p in the represen- 
tation n^, 

(z + j — a + a a + \ 

ro- + -^ ■_,+-■ (3 ' 34) 

(2j - a+a ) a z - j + a+a I 

Its matrix elements depend on two parameters z and j. It is convenient to define new 
variables 

z ± = z±{ 3 + \) . (3.35) 

Consider the limit 

z — )■ <x>, j — > oo, -2+ = fixed , (3.36) 

and define 

L +<-> = & (J A) ^+ - i - \) = (^ f ) • (3-37) 

where h is defined in (I3.26p . Similarly, for fixed z_ define 

£.(«.)= ]JmL,(*. + J + i)(f ^M ^; h ). (3.38) 



The particular transformations of the form (13. 33 p . used in the definitions (I3.37P and ([3T 



do not affect the validity of the Yang-Baxter equation (13. ip . Therefore, the new operators 
L±(z) will automatically satisfy (13. ip in virtue of the commutation relations of the oscillator 
algebra. Of course this may also be verified by direct, elementary calculations. It should be 
mentioned that these two L-operators are not really new objects, and appeared earlier in 
different contexts. They are e.g. known to yield the Lax operators of the so-called "discrete 
self-trapping chain" [251 EZ] • 

A typical problem frequently arising in the theory of integrable system is the "fusion" 
of different solutions of the Yang-Baxter equation. This is a standard way to understand 
relationships between various solutions and to obtain new ones. It turns out that in our 
case the consideration of the product of L + (z) and L^(z') does not lead to any new non- 
trivial solutions. It however allows to discover remarkable factorization properties of the 
operator Lj(z), which will be explained below. 

Consider the direct product H® 2 of two oscillator algebras (I3.26p . Let {af,hi} and 
{a 2 , h.2} denote the two associated mutually commuting sets of generators. Below we will 
use the following similarity transformation acting in this direct product 

5± = e a+a 2 " a ± e -a+a 2 " ^ k = 1,2. (3.39) 
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Obviously, the transformed operators af 2 obey exactly the same commutation relations as 
af 2 . In particular, the operators a^ 1 commute with a^. Explicitly one has, 

af = aj~ — aj , a^ = a^ , (3.40) 

a 2 = a 2 i a 2 = a i + a 2 • (3.41 J 

Consider the product 

^^tf'w^w^' ,_ a + v,)( 2 -;" 2 ?)• ^ 

where the superscripts (1) or (2) indicate that the corresponding L-operators belong re- 
spectively to the first or the second oscillator algebra. By elementary calculations one can 
bring this product to the form 



1 a+ \fz + -h 2 a+\ +aj -/l 0\ (z + 3) 3J 

ar 2 - + h liy H -a 2 i) ' V a r i ; V J t z ~ J " 



(3.43) 

where h 12 = a^ 2 a^ 2 + | and the sl(2) generators J|, J+ and J~ are realized as in f )3.3ip 
but employing the operators a^, 

J" = a+ , J+ = (2 j - a+ aj) aj , J 3 = j - a+ a^ . (3.44) 

Introducing operator valued matrices 

*-(-)< M-)' (345) 

one can re-write (13.431) in a compact form 



C(z) = L W ( z -j-l)L ( ?(z+j + ±) = e<^ B W hf (z) e~ a ^ , (3.46) 

where the superscripts (1) and (2) have the same meaning as in f)3.42p . Similarly one 
obtains 

X(z) = L { l\z + j + \) L<?(z - j - \) = e a ^ L?\z) B® e -** , (3.47) 

Writing the last identity in full, one gets 

where the sl(2) algebra is now realized as in (I3.32p with a 1 

3 J = a+ (2 j - a+ a7) , J+ = a7 , 3) = j - a+ a7 . (3.49) 

15 



3.4 Construction of the Q-operators 

We are now ready to explicily define the operators Q±(z) as transfer matrices. To do 
this we will use our general definition (13. 8p of a transfer matrix for the L-operators L±(z) 
defined in (13.371) and (13.381) . Solving (13.91) one gets the same boundary operator in both 
cases, 

D± = e~ i(ph , (3.50) 

where h is defined in (13.261) . In this way we define operators 

Q±(z) = Z- 1 e^ z Tr T (M±(z)), (3.51) 

where 

M±(z) = e~^ h L±(z) <g> L{z)± ® ■ ■ ■ ® L±{z) . (3.52) 

S v ' 

L— times 

Note that we have changed the normalization of (I3.5ip in comparison with ( 13. 8 \ by intro- 
ducing z-dependent exponents and a constant factor 

Z(0) = Tr^(e^ h )=-^ T . (3.53) 

1% sm | 

where the subscript J 7 indicates that the traces are taken over Fock space (13.271) . By 
construction the operators (I3.5ip will automatically commute with T(z) as a consequence 
of the Yang-Baxter equation ( 13. ip and properties (13 .4p and ( I3.9p . 

[T(z),Q ± (z')} = 0. (3.54) 

The operators Q±(z) can be obtained from each other by negation of the twist 0. It is 
not difficult to show that 

Q + (z,(j)) = 1Z Q_(z, — 4>) 71, 71 = a x <8> a x <g> • • • <8> a x , (3.55) 

where 7L is the spin reversal operator in the quantum space. 

Let us now prove the factorization equation (13.181) . Consider the quantity defined by 
the first equality in ( I3.46p . 

C{z) = L W {z- j -\)Lf{z + j + \) . (3.56) 

where superscripts (1) and (2) have the same meaning as in (I3.42p . This is a product 
of two L-operators with different auxiliary spaces (each being a copy of the Fock space), 
which, of course, satisfies the Yang-Baxter equation itself. It can be regarded as a single 
L-operator, whose auxiliary space is the tensor product of two Fock spaces. Therefore, one 
can apply our general formula ( 13. 8 j) to define the transfer matrix 



T 



c (z) = Tr TxT% (e-'*( hl+h2 > £(z) ® C(z) 8-8 C(z\ ) , (3.57) 



V 

L— times 
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where we have substituted 3± defined in (13.501) . It is not difficult to see that using (I3.56P 
one can rearrange factors under the trace such that (I3.57P reduces to a product of two 
operators (I3.5ip . 

T c (z) = Z{<pf e-^ + b Q_( Z -j -I) Q + (z +j + l). (3.58) 

Similarly, using the second expression for C(z) from (I3.46P one obtains 

1 



T £ (z) = e -^(J+5) J B _ Tf(z), (3.59) 

J B _ = Tr T (e-^ h B_ <g> B_ ® • ■ ■ ® B_ ) . (3.60) 



'3 

where 



L— times 



Note that the presence of the similarity transformation in (I3.46P does not affect the calcu- 
lation of the trace since the operator e ai &2 commutes with hx + h 2 in the boundary twist. 
The matrix B_, defined in ( I3.45p . is triangular. The calculation of the trace in (I3.60p is 
trivial, since _B_ depends only on a - , it gives 

T B _ = Z(<t>)I, (3.61) 

where I is the unit operator in the quantum space. Combining everything, one arrives at 
the factorization formula (I3.18p . though with swapped order of Q + and Q_ in the product. 
Repeating the same reasonings, but this time starting from the equation (13.471) . one obtains 
(I3.18P exactly as written, which proves the commutativity of Q + and Q_, stated in (13.191) . 



Let us stress that one should explicitly show that also the last equation in (13.191) is 
satisfied by the operators constructed above. This can be proven in the usual way starting 
from the Yang-Baxter equation in C 2 <g) osc £g> osc. The existence of an intertwiner in this 
case has been shown in [27] in the course of a study of the DST chain, which is a certain 
bosonic hopping model. The further investigation of this and related issues will be reported 
elsewhere [28] . 

For the benefit of the reader, we will present for small chain lengths L = 1, 2 the explicit 
forms of the finite twist operators T, Q± as well as their eigenvalues in appendix [Bj 

3.5 Bethe equations without a Bethe Ansatz 

It is interesting that our factorization formulas (13.431) and (I3.48P allow to solve the twisted 
Heisenberg model without the somewhat tedious Bethe Ansatz technique, be it a coordi- 
nate, algebraic or functional Bethe "Ansatz"! In fact, no Ansatz (German for "try and see 
whether it works") for some wavefunction is ever made. One instead derives the funda- 
mental operator relation (I3.20p as done in the last section. The full hierarchy of operatorial 
Baxter and fusion equations immediately follows, as we explained in the course of the dis- 
cussion of (I3.24p . In particular, one immediately derives the operatorial Baxter equation 
(I2.32p from (I3.20p with auxiliary spin j = \- 
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As shown above all T- and Q-operators belong to the same commuting family and 
therefore can be simultaneously diagonalized by a ^-independent similarity transformation. 
Thus the eigenvalues will have the same analytic properties in the variable z as the matrix 
elements of the corresponding operators. By construction all T-operators are polynomials 
in z, while the Q-operators are polynomials multiplied by the simple ^-dependent exponents 
(see (13.511) ) so that their eigenvalues are exactly of the form (12.211) . (12.171) . Furthermore, 
all the functional equations can be considered in the basis where all operators are diagonal 
and can be replaced by their eigenvalues corresponding to the same eigenstate. In other 
words the functional equations can be treated as scalar equations for the eigenvalues. For 
instance, substituting (I2.2ip . (12.171) into (I3.22p one obtains 

e^A+(z + §) A_{z - |) - e-&A.(z + \) A + (z - \) = /(</») z L . (3.62) 

Let Zk be a zero of A + (z), i.e, A + (z k ) = 0. Using this fact in the last equation one obtains, 

e~**A4z k ) A + (z k - 1) = (z k - |) L , e&A+(z k + 1) A.(z k ) = (z k + ±) L , (3.63) 

Dividing these equations by each other one arrives at one set of the Bethe equations (I2.18p . 
The other set is derived similarly. 

We would like to mention one simple but important corollary of the quantum Wronskian 
relation (I3.62p . It concerns exact strings, i.e., groups of zeroes {zi, Z2, ■ ■ ■ , zi} equidistantly 
spaced with the interval 1, 

z k+l = z k + l, k = l,...,£-l, (3.64) 

where £ is called the length of the string. It is obvious from (I3.62p that neither of the 
polynomials A±(z) could have exact strings of the length greater then two. Indeed, the 
only such zeroes can be at z = ±|, thus A±(z) might only have (possibly multiple) 2-strings 
at z — ±|. The same analysis extends to the zero- field case as well. Note, in particular, 
that the long strings on the imaginary axis of the variable u = iz, shown in fig. [TJin section 
14.51 are not exact. 



4 Removing the Twist: The — >- Limit 

In the previous chapter [3] we have used a regulator <p to obtain finite quantities when 
taking traces over infinite dimensional oscillator spaces. It is interesting that this natural 
quantity corresponds precisely to the spin chain with magnetic flux discussed in the 
review chapter HJ In the present chapter we will deal with the operation of taking this 
regulating flux away. This involves a significant symmetry enhancement to a global SU(2) 
symmetry of the spectrum of the Heisenberg chain, "weakly" broken by the twist. Moreover 
we would like to make contact with the results of [18] . see also the attempt in [29], as well as 
the findings in [23]. It is interesting to explicitly work out the operators Q± for small spin 
chain lengths L. This is done in appendix [B] In line with their construction, the operators 
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Q-t will diverge as the regulating is removed, and we will need to "renormalize" them. 
The reason is simply that the sum over oscillator states ceases to converge, as maybe 
explicitly verified on the examples in appendix [B] Let us begin by studying the structure 
of divergencies of T^ and Q±. 

4.1 Divergencies in the zero twist limit 

It is interesting to understand the singularity structure 15 ! of T, T + and Q^ in the limit 
where the twist <fi is sent to zero, and su(2) invariance is recovered. In preparation, we 
will prove in appendix [C] that Tj(z\ <p) is finite in the zero- twist limit for arbitrary complex 
values of z and j, which is a priori not obvious. It is trivially finite for 2j e Z + , since in this 
case there is near complete cancellation between two infinite sums, the remainder being 
a sum over 2 j + 1 terms. As for T^(z; 0), one can immediately say that its eigenvalues 
diverge 16 ! at most as <p~ L in the zero-twist limit. However, we found evidence that the 
actual divergence is much milder, namely 17 ! 

T+(z;0)^<T 2 ^,j), (4-2) 

with some matrix 0(z,j). Here J is defined as 

J(J + l)=$ ot . (4.3) 

This being a quadratic equation, there are two solutions for J~, and one has to choose the 
one such that the spectrum of 1J belongs to Z> . 

In any case, we already showed that the operator T~-"(z) is a "composite object", and 
thus perhaps not too interesting. Much more relevant is the zero-twist limit of its fun- 
damental constituents Q±(z±). As we discussed in the review chapter |2J we do expect 
these operators to become singular, as Bethe roots at infinity physically indicate su(2) 
descendent states. Their naive divergent behavior is <p~^ ±s \ but their actual divergence 
is again empirically found to be much milder 

Q±(z±;4>)^r (J±S3) OHz±)- (4-4) 

In particular their action on lowest and highest weight states is finite in the zero-twist limit 
for Q + and Q , respectively. One can check on small-length examples that ± (z) defined 



15 We do not have proofs for the singular behavior of the operators T + and Q* discussed in the following. 
Our statements are based on observations on spin chains of small length L = 2, . . . , 5. 

16 Note that 

OC 

Y,e~ Xn n l -{l-^-x~ {l+1) (4-1) 

n=0 

in the x — > limit. 

17 It is not clear to us whether 0(z,j) can be meaningfully interpreted as a "renormalized" transfer 
matrix. E.g. we do not know whether this operator is su(2) invariant. An alternative, more symmetric 
definition would be Tt(z; cf>) — > (f)~ J 0(z,j) <f)~ J . 
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in this way is not su(2) invariant. An su(2) invariant operator Q(z) may be denned as 

Q(z) = lim iJ ^f )] (T<P)^ J±s3) Q±(z-, <P)(T<P)^ J±S3) , (4.5) 

where it may be observed that the eigenvalues of Q are polynomials in z normalized as 
z 2J + . . . . We see that both Q± lead to the same Baxter operator Q, see also [23]. This is 
technically easy to understand since Q± are related to each other by a spin-flip combined 
with reversion of the twist ~* ~ $>■ An su(2) invariant operator like Q at = is 
obviously invariant under this transformation. However, this means that in the limit we 
are "loosing" one of the two linearly independent operators of the twisted case. 

4.2 Q-operator in the zero twist limit 

In this section we will propose another way to define the Q-operator in the zero-twist limit. 
This second way is certainly less explicit than (14. 5p . but is more suitable for rigorously 
proving finiteness, as well as su(2)-invariance of the resulting Q-operator. The renormalized 
Q(z) can be very naturally defined by "renormalizing" Q,±(z), namely 18 ! 

Q zo (z) = lim Q ± (z; 0) Q^; 0). (4.6) 

The resulting Q still carries a dependence on the "subtraction point" z . In this sense (14. 6p 
defines a one parameter family of Q-operators. Q Zo (z) in (14.61) satisfies Baxter's equation. 
This is easily seen since Baxter's equation is linear in Q, and, at fixed values of 0, Tj(z\ 0) 
and Q±(.z;0) commute among themselves. It is important to stress that the Q-operators 
defined in (14.61) are su(2) invariant. The proof of this fact is left to appendix [D] su(2) 
invariance is the underlying reason why in ( I4.6P both Q + and Q_ lead to the same Q. Let 
us also point out that the obtained Q-operator will turn into the identity operators for 
z = zq, an observation which will be used later. 

The fact that (14.61) is a finite quantity is easy to understand. Indeed one can safely 
discuss the finiteness issue at the level of eigenvalues. The similarity transform that diago- 
nalizes Q± is well defined for any (including = 0), being the same similarity transform 
that diagonalizes Tj(z] 0). Any given (fixed) eigenvalue of Q±{z) will be of the form 

g ± (z,0) = 0- fc (g ± (^) + O(0)), (4.7) 

for some integer k. The ratio 

q±{z ' A) (4.8) 



<?±(^o;0) 



18 Notice that Q 2o (z) = Q(zq)~ 1 CI(z). Something similar appears in Baxter's original work [7]. In [7] 
Baxter built two solutions to his equation, termed Ql,Q_r. These two matrices do not commute with 
the T matrix and do not commute among themselves (they satisfy Ql{u)Qr(v) = Ql{v)Qr(u)). The 
operator Q(u) = Qfl(w)Q i? 1 (uo) = Ql(-u)Q^ 1 (u ) is introduced to cure these two problems. We find this 
similarity quite interesting. 
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will be finite in the limit. Of course (14. 6p requires det Q±(z ; (f>) ^ 0. To be precise, z 
should not be a Bethe root for all values of in a small open domain in the vicinity of 
= 0. The zero twist limit of (I4.8P is a polynomial in z normalized to be one for z = zq. 
The degree of this polynomial can be worked out using su(2) invariance of Q Zo (z) and the 
fact that Q + and Q_ are finite (because the Bethe roots stay finite) in the zero twist limit 
for, respectively, lowest and highest weight states. Moreover notice that the exponential 
in (I2.2ip can be simply replaced by one in the limit (14.61) . 

We finish this section by noticing that any two Q zo (z) defined in (14. 6 p are linearly 
dependent since 

Q Z0 (zi) Q z >(z 2 ) - Q Z0 (z 2 ) Q z , o ( Zl ) = . (4.9) 

In the next section we will discuss how to obtain another one-parameter family of solutions 
to Baxter's equation, which will turn out to be linearly independent to the ones we just 
obtained in (14. 6p . 

4.3 P-operator (as a T-operator) in the zero twist limit. 

In the previous section we identified a one-parameter family Q zo (z) °f solutions to the 
Baxter equation in the zero twist limit. As already discussed in section [2], Baxter's equation 
possesses two linearly independent solutions. Accordingly, let us now obtain a second, 
linearly independent family of solutions, which we will denote by P Zo (z). While constructed 
in a different way, it is identical to the one obtained in [29J, and is based on the same 
observation 19 !. We will relate this second one-parameter family of operators to Q±(z;0). 
The observation made in [29] is that 

PJ^eT^-i(^), (4.10) 

2 

satisfies Baxter's equation. z is then the label of the one-parameter family. Notice that 
in this equation there is no twist and P Zo (z) is an su(2) invariant operator. 

For P Zo (z) to be defined for any complex z, the T-operator has to be constructed 
for any complex value of the spin. Different approaches to this problem are reviewed in 
appendix O They include the construction proposed in [29], the trace functional approach 
introduced in [3U] and the approach of the present paper (see also [25]). 

It is instructive to explicitly write (I4.10p using the fundamental relation (13.201) 

P Z0 (z) = lim r\<f>) (Q_(*; 0) Q + (z ; <P) ~ Q-(*b; <P) Q+(z; 0)) , (4.H) 

and to take a closer look at the mechanism by which the eigenvalues of P ZQ {z) become 
polynomials of the degrees found in [18]. A proof is encoded in the zero twist limit of the 
generalized Wronskian equation analyzed in the next section, here we will merely give some 
heuristic argument. In contradistinction to the Q zo (z)-operator, the P 20 (z)-operator has 



19 To avoid confusion, we note that our P was called Q in jjJS]. However, its eigenvalues coincide, up to 
overall state-dependent normalization, with the polynomials P in [18j . 

21 



been denned as the difference of two quantities. Finiteness of ( 14.111) is then a consequence 
of the cancellation of many divergent contributions from the two terms. This means that 
in the expansion of ( 14. lip the exponential factors appearing in (I3.5ip and ( I2.2ip play 
an important role. Their major role is to increase the degree of polynomiality of the 
P-operator as compared to the naively expected one. 

As opposed to the one-parameter family of Q-operators, any two P-operators at dif- 
ferent values of the parameter z are linearly independent. This is the underlying reason 
behind the appearence of an extra root in the P-operator's eigenvalues: in the zero twist 
limit the P-operator contains the Q-operator but not vice versa. Let us next analyze this 
in more detail. 

4.4 Generalized Wronskian for P and Q 

Thanks to the insights gained in the previous two sections, we will now study the relation 
between the Q-operators and P-operators given in ( 14 .6p and ( 14. lip . It is interesting to 
notice that they can be consistently defined for any value of the twist 0. Inspired by (14. 6p . 
( 14"TT|) . let us define 

(4.12) 

(4.13) 

The explicit form of /(0) is not important here, the important thing is that it is the same 

function that appears in the fundamental relation. As already shown, both Q ZQ (z; 0) and 
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Pz (z',4>) satisfy Baxter's equation_| and are finite in the zero twist limit. Here we want 
to stress that they satisfy generalized Wronskian relations. It can be shown by direct 
substitution, using the fundamental relation (I3.20p . that 

T J (z;^ = P Z0 (z + j + l; ( p)Q Z0 (z-J-^(p)-P Z0 (z-j~^ ( p)Q Z0 (z + j + ^^. (4.14) 

A question that naturally arises is the meaning of the parameter zq. It is interesting to 
notice that Tj(z] 0) in (14. 14p does not depend on the choice of Zq. On the other hand, in 
the very same way as (I3.23p . (I3.20p have an intrinsic sl(2) invarianca 21 ! which rotates Q± as 
a doublet, Tj(z; 0) in (14.141) is left unchanged under rotations^ 2 ] of Q,z ( z ; 0) and P Zo (z; 0). 
It turns out that these two invariances of (14.141) are not unrelated. It is easy to see that 



P-iteM __ (&$\ (Pzo(z;<t>)\ „ (QMA) Pz' (z -A)\ [P zo (z 



Qz' (z;(f))J \*f 6 J \Qz (z; 0) J ~ V N ^(^o;0) Qz> (z ; 4>) J \Qz {z; 



(4.15) 



20 Baxter's equation is the same for all these operators and does not contain the twist explicitly. Its 
information is hidden in the analytic structure of the solution we want to obtain. 

21 Note that this sl(2) is not identical to the manifest su(2) of our Heisenberg magnet under study. 

22 In appendix El we review how Q 2o (z;0) and P 2o (z;0) transform under s((2) rotations of Q± as a 
doublet. The fact that a u(l) £S> Z 2 inside fil(2) acts freely on Q Zo (z; (j>) and P Zo {z; <f)) results in equations 
(|4.12|) and (|4.13p not being invertible for Q ± . 
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where 

N 4 (Z O ;0) = — (q^.^Q^.^) " Q_( 2b ;0)Q + ( Z /;0)J " (416) 

We have a 5 — (3 7 = 1, which is consistent with the sl(2) rotation. This identity is a special 
case of the following interesting relation 

N 20 (4 <f>)P z ,(z; <P) = Q Z0 (z; <P) Q z ,(z'; 0) - Q Z{) (z'; 0) Q 4 (z; 0) . (4.17) 

Note that N Zo (z ; ;4>) vanishes in the zero twist limit, in line with equation (14. 9p . Let us 
stress the implications of this fact from the point of view of equation (j4.15p . In the zero 
twist limit varying zq corresponds to triangular rotations of P and Q. In particular, it 
means that the Zq freedom corresponds to the normalization of Q, but does not affect its 
Bethe roots. For P the freedom in the parameter zq corresponds to adding the Q-operator 
to the P-operator. From the point of view of the Bethe roots distribution it corresponds 
to the position of one of the roots of P. It is indeed obvious from equation (14.121) that 
P 20 (z ; 4>) = 0. The positions of the other Bethe roots then depend parametrically on the 
position of the root z . Note that the eigenvalues of (14.121) and (14. 13j) for the highest weight 
state in the zero twist limit will be 



M M L-M 

z — Zk 



Qz {z) = TT k - , P Z0 (z) = (2M + 1) TT (z - z k ) ]"[ (z - z k (zo)), (4.18) 

Zn — Zh 
fc=l u fc=l fe=0 

where zq(zo) = zq. The eigenvalues of all the other states follow from su(2) invariance. 



4.5 Numerical results on the root distributions of Q(u) and P(u) 



In the previous sections^ we obtained for twist (p = a one-parameter family of a pair of 
operators {Q(w), P(w)}. On the level of eigenvalues, this free complex parameter changes 
the overall normalization of the polynomial eigenvalues Q(u) in ( I2.30p . and both the nor- 
malization and the position of the roots of the polynomial eigenvalues P(u) in (12. 30 p . Let 
us numerically study the positions of the roots of these two complementary polynomials 
for various interesting classes of states. Here clearly the overall normalization does not 
matter, and we can use the parameter to shift one of the L — M + 1 roots of P(u) to any 
position in the complex u plane we fancy. 

We will focus here only on the case where both L and M are even. Let us start 
with the polynomial Q(u). From the definition it is a polynomial of degree M with roots 
being the solutions of Bethe equations (12. 7p . With every Bethe root u^ we identify an 
integer n^, the so-called mode number, which labels different solutions. Every solution is 
then characterized by its set of mode numbers. To keep the discussion as transparent as 
possible we will focus only on specific choices of this set. Namely, we will consider two-cut 



23 Recall that for the numerical studies of this section we will return from the z-plane to the Wick-rotated 
u-plane via (|2.3ip . 
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Figure 1: The two-cut root distribution of the polynomials Q(u) (purple squares) and P(u) 
(blue dots). Here L = 60, M = 22 and \n k \ = 1 for all k. 

solutions by taking symmetric distribution of roots on the complex plane with all mode 
numbers equal \n k \ = n for k = 1, . . . , M. The method of finding such solutions can be 
found in [31] (see Appendix|F]for details). With the algorithm given there we are however 
only able to find solutions of ( 12. 7p when restricting to H- < |, i.e. below half-filling. In 
order to then also find the dual roots from (12. 9p . which corresponds to the polynomial 
P(u), we have to find a different method. The idea is to use the generalized Wronskian 
relation (j4.14p with = and j = 0. Then the left hand side of (j4.14p is exactly u L , and 
the problem of finding the polynomial P(u) reduces to a linear problem for its L — M + 2 
coefficients. Having the polynomial, finding its roots is a trivial numerical task. 

Using the methods presented above, we were able to find many different configurations. 
One typical example is given in figure [U see also [32] and [33]. Here the length is L = 60 and 
the number of Q-roots is M = 22. As we mentioned before, there is a freedom in choosing 
the position of one of the roots of the polynomial P(u). We took it to be uq = Q 24 l . As was 
pointed out in [31], the positions of the roots of the polynomial Q{u) are well described by 
the predictions from the thermodynamic limit - they lie almost on a "spectral curve" . On 
the other hand, the roots of the polynomial P(u) are located on a dual configuration of 
cuts (see fig. 40 in [3T]). It means that Q and P establish two different ways of slicing up 
the complex plane, with cuts originating from the branch points given by thermodynamic 
calculations. 

We will end our discussion of Q- and P-roots with a few comments on the mode numbers 
rik relevant to the above examples. In order to define these, the untwisted Bethe equations 
( 12. 7p need to be rewritten in logarithmic form 

Llog 2 - + 2mn k = V log — 3 - :, 4.19 

u k - \ j^ k u k - Uj -i 



24 In the Appendix IF1 we will also show results for other choices. 
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where the summation is over all roots of either the Q- or the P-polynomial. We can 
immediately check that for the Q-roots presented in fig. [T]the mode numbers are +1 for 
the roots with positive real part and —1 for the ones with negative real part. The same is 
true for the P-roots, but in this case there are also other roots which are purely imaginary 
and form a so-called condensate, which is nothing but a giant string in the standard 
terminology of the Bethe Ansatz (to be precise, in the case of fig. [T] it is a combination of 
two strings shifted by |). As explained at the end of section 1331 the strings are not exact, 
however their roots very closely approach the (half) integer positons on the imaginary axis 
(the roots at u = ±| are exact). The mode numbers for the roots in the strings, are 
not well defined due to the ambiguity of defining logarithms of real negative arguments, 
appearing in r.h.s. of f !4.19|) . 

5 Relation to Earlier Work 

5.1 Historical note. 

The concept of the Q-operator had been introduced by Baxter in his seminal paper [7] on 
the symmetric 8-vertex model and further developed in this context in [M| [35] and more 
recently in [SUEZ]. Baxter's idea received a major boost in jB], in an attempt to understand 
the integrable structure of conformal field theory, which is related to the continuum limit 
of the integrable XXZ model. Renewed attention to the Q-operator arose independently 
in the context of the discovery of non-compact integrable Heisenberg magnets in the high 
energy scattering problem of QCD [38, [39] . The emergence of integrability in the spectral 
problem of planar M = 4 and in AdS/CFT [TJ [2J El H], not yet understood beyond the 
one-loop level, is the motivation behind our fresh look at the Q-operator in this article. 

In the present chapter we will very briefly review what was previously known about 
the Q-operator for the XXX spin j Heisenberg magnet in order to compare with our 
findings above. In the comparison between different approaches to the problem it seems 
to be important to distinguish between compact magnets (2j e Z> ) and non-compact 
magnets, where typically either 2j = —1, —2, . . . (discrete series representations of su(l, 1) 
or sl(2, R)), or else j takes on certain continuous values (principal and supplementary 
series representations). The quantum space is indeed very different in the compact and 
non-compact cases, as it is finite dimensional in the former, and infinite dimensional in the 
latter for any length L of the chain. 

5.2 Non-compact magnets 

In the case of non-compact magnets two main ways to built Baxter's Q-operator(s) have 
been proposed. These are related, but this fact has not yet been properly investigated. 
The first way is based on the connection between Backlund transformations in the theory 
of classical integrable systems and Baxter's Q-operator for quantum integrable systems. 
Sklyanin et al. (see jlQ] and references therein), inspired by [H] , proposed a way to build 
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the Q-operator as an integral operator for models governed by the s [(2) invariant R- 
matrix. It is constructed as the trace over some monodromy built from Lax operator 
intertwining operators connected with the Toda, DST (discrete self-trapping) and XXX 
models, respectively. An explicit construction for the XXX model has been proposed in 



The second way is based on some factorization properties of matrix elements (related 
but different to operator factorization of this paper) of i?-matrices and L-operators, which, 
apparently, was first employed in (33] in the context of the chiral Potts model [H] SSI US] • 
This approach allows to explicitly calculate matrix elements of the Q-operator as a simple 
product involving only two-spin factors. The same method was later used in J5T] for the 
Toda model then in a series of papers [25j IS] studying various aspects of the factorization 
of i?-matrices, associated with the XXX model. It would be interesting to clarify the exact 
relationship (if any) between the factorization in these works and the one used in this 
paper. 

5.3 Compact magnets 

Among the compact magnets XXXi/ 2 is the most studied, beginning with Bethe's original 
work [13J. Here we will restrict to this case. While a lot is known about the Q-functions, 
i.e. the eigenvalues of the Q-operator (s), much less is known about the construction of the 
latter as a trace over a suitable monodromy matrix for XXXi/2- 

In [T5], inspired by [9], Baxter's equation, at the level of eigenvalues, had been consid- 
ered. For a given eigenvalue of the T-operator acting on some h.w.s, two linearly inde- 
pendent polynomial solutions of Baxter's equation were found. These two solutions were 
called P and Q in [18], and are generically distinguished by their degree of polynomiality in 
the spectral parameter. It was furthermore shown that the lower degree polynomial Q can 
always be added to the higher degree polynomial P, leading to a one parameter freedom 
in the root distribution of the higher degree polynomial. 

Subsequently Pronko built a one-parameter family of Q-operators in [29], which we 
denoted in section fl~3l P-operators. However, as explicitly written in equation (I4.10p . the 
one-parameter family of P-operators is nothing but the set of all T-operators after a linear 
redefinition of their spectral parameter and spin. Thus the construction in [29], following 
[7], finally boils down to the construction of Tj(u) for general, complex values of the spin 
j G C. As discussed in appendix O there are various seemingly different but ultimately 
equivalent procedures to achieve this. 

In [1H] and [27] Pronko tried to build two "basic" Q 1 * 1 instead of a one parameter family 
of Q-operators. He succeeded in doing this for the Toda chain and the DST "discrete self 
trapping" chain, respectively. Still it is important to keep in mind that his construction, 
apart from the fact that it was applied to a different class of models, remains conceptually 
different from the one proposed in this work. The idea in [291 l4~8l [27] is to build some oper- 
ator satisfying Baxter's equation. In contradistinction, in the present paper we aim at the 
fundamental relation (13.201) . The factorization formulas (13.431) . (13.481) are the crucial step 
in doing this. The relation between the two constructions deserves further investigations. 
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6 Conclusions and Open Problems 

In this work we have demonstrated that the complete understanding of the "simplest" of 
all integrable quantum models, the compact spin-| XXX Heisenberg model, requires the 
introduction of non-compact oscillator representations. After this is done, a very inter- 
esting factorization of the compact quantum Lax-operator takes place, cf. (I3.43J) . (j3.48p . 
This fact may then be used to quickly find a very elegant solution of the model in terms of 
the full set of operatorial (as opposed to mere equations for their eigenvalues) functional 
relations between a large class of commuting families of transfer matrices, which natu- 
rally includes two linearly independent Baxter Q-operators. An added benefit is that all 
transfer matrices, including the Baxter Q-operators, are explicitly constructed as traces 
over suitable monodromy matrices. The latter is particularly important, as the analytic 
structure of the solutions to the hierarchy does not have to be guessed, but is clear from 
the construction. In this sense we consider our approach a first sketch of what needs to 
be done in order to understand the solution of the AdS/CFT spectral problem proposed 
in [5] , where the analytic structure of the solutions of the infinite set of functional integral 
equations obtained from the TBA approach has to be guessed and then imposed by hand. 

Apart from this technical advantage of a direct construction on the operatorial level, we 
believe that there will eventually be an interesting physical interpretation of our procedure. 
Baxter Q-operators enter the T-system of AdS/CFT as "boundary values" . Furthermore, 
it is clear that the auxiliary channel in which the bosonic excitations we need for our con- 
struction are propagating is closely related to the "mirror direction" in the TBA approach. 
A similar analogy was recently pointed out in [49]. We believe that these kinds of funda- 
mental short-range excitations are suitable for replacing the long-range "mirror magnons" 
of AdS/CFT. 

An important feature of our approach is that a magnetic flux (horizontal field, or twist) 
has to be applied in order to ensure convergence of the Baxter Q-operators. This suggests 
that in some sense the theory with the flux is more natural, and the untwisted case is a 
somewhat singular limit. It is possible that this feature also lifts to the AdS/CFT case. 
As was stated in the introduction, in AdS/CFT the twisted theory also appears to be 
integrable [19] , but much less is known about it. In particular, its finite-size structure is 
very puzzling, see e.g. [20] ■ In any case, from our point of view it is still not clear whether 
one can construct Baxter Q-operators in the absence of horizontal fields directly as a trace 
over some monodromy matrix. 

Clearly our approach should be extended to the compact si(n) and the supersymmetric 
sl(n\m) cases, as well as the non-compact magnets, with the goal of treating the symmetry 
algebra su(2, 2|4) of AdS/CFT [23]. In particular it will be interesting to relate our method 
to the Baxter Q-operator construction developed for non-compact spin chains in [23| |4"T]. 

Finally we would like to mention yet another interesting research direction in the theory 
of Q-operators. It concerns their connection to the spectral theory of differential equations 
[51~| [52], [53]. It recently found new applications in the theory of BPS states of a large class 
of d — 4, J\f — 2 field theories [54] , as well as in the (mysteriously) closely related theory of 
classical string solutions for certain strongly coupled Wilson loops in M = 4 gauge theory 
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J. Lastly this spectral theory also appeared in the description of the massive Sine(h)- 
Gordon model [55] • For the XXX model the connection to differential equations arises in 
the c = 1 CFT limit [57], which is closely related to the Kondo model |58j . 
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A Functional relations 

Using the defining relations for T and T in f !3.23|) and (13. 20 p . one finds 

TAB = f(4>)TA-B-l(±t£). (A.l) 

2 

A more explicit form of the left equation in (I3.24p reads 

Tcb Qa = Tab Qc + T CA Q| . (A.2) 

Choosing 

A = z, B = z + z'+j + \, C = z + z'-j-\, (A.3) 

one gets 

Tj{z + z>) Q ± (z) = Ti & _ y _i } (z + l(z'-j- A)) Q ± (z + z' + j + §) + 

Ti ( . + ^_i } (z + \{z' + j + I)) Q ± (z + z'-j-l). (A.4) 

Analogously the right equation in f)3.24[) reads 

TcbTad = Tab Tcd + Tca Tbd ■ (A. 5) 

With the choice 

A = z' + j' + ^ D = z'-j'-l, C = z + 3 + \ B = z-j-\, (A.6) 
one gets 

T J (z)T,(,Q = T ^ + ,_ / _ 1 ( Ww + i ) t ,,,^,, ( *+*-i-i'-A + 

2 v ' 2 V ' 

T *^±f±i (^ ± ¥= i ) T ^±2±l i*- 1 ^-) ■ (A.7) 
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B Explicit results for small chain lengths 



In this appendix, where we revert to the notation 25 ! with the spectral parameter u = iz, 
cf. (12.31 j) . we will present the results for Q_(it), Q+(w) and T(w) for the smallest possible 
chain lengths^ L = 1,2. All matrix elements of Q_(w), Q+(w) are easily computed 
from f)3.5ip . For the convenience of the reader we will present the "Wick-rotated" (i.e. 
written in the w-plane) Lax operators for the Baxter Q-operators. Recall also the definition 
z± = z ± (j + |) in (I3.35p . They read, with h = a + a" + \ as defined in (13.261) . 

L"(w)= ( l * + .. > ) . (B.l) 

and 

L + H=( M _" a - h *f ) ■ (B.2) 

The definitions of L given above differ from the definitions L± given in (I3.37p . (I3.38P in 
the main text. The two definitions are related in this way 

IT(«) = (j ° ) 1*4*) , L+(^) = (o i ) L + {z) . (B.3) 

The conventions used in this appendix are such that the eigenvalues of Q± have the form 
exp (±|m)(m p + . . . ), where p is some integer. Notice that the position of the diagonal 
matrix (left or right) in (IB.3|) can be compensated by an oscillator algebra automorphism. 
In particular this means that left and right multiplication became equivalent after the trace 
is taken. This in nothing but restating the U(l) invariance of Q±. The Baxter operators 
(I3.5ip become in the w-notation 

Q ± ( M; <P) = T J ( ^ fe) T^ (e-^Lt (u) (g> • • • (g) Lf (ti)) . (B.4) 

The matrix elements of the transfer matrix T(w) are then obtained, for j = ~, from (I3.20p . 
It is easy to check that the direct construction based on (12.241) . (I2.25P leads to the same 
results. We will also write down the L — 1, 2 eigenvalues Q-(u), Q+(u) and T(u) of these 
operators. 

It is easy to check that the operators below indeed satisfy the j = and j = \ cases of 
the Wronskian relation (I3.20p translated to u-space 

2i sin ^u L = Q + (u + l -) Q_( M -'-)- Q + (u - % -) Q_( M + % -), (B.5) 

2isin|T(M) = Q+(u + i)Q_(M-i) - Q + (u -i) Q_(w + i). (B.6) 



25 The reason is that the Bethe equations are most commonly written in the u-plane. 
26 Note that while it is not possible to define a nearest neighbor Hamiltonian for a L — 1 "chain" , the 
definitions (|2.24j) . (|2.25p for the transfer matrix and p.51[) for the Baxter operators still make sense. 
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C Analytic continuation of the trace 

In this appendix we will review three ways of defining traces for complex spin. The first 
way, used in our construction in the main body of the paper, involves introducing a twist 
as a regulator and then takes (f) — > while subtracting infinities. The second way has 
been introduced in [30] under the name "trace functional". The third way has been used 
by Pronko in [29]. We will now show the equivalence of these three constructions. In the 
case of the first two, this had already been noticed in [23J . 

Let us first fix the set-up for the analysis. Denote by O any one of the 2 L x 2 L matrix 
elements of the monodromy matrix, see (I3.12p . O is some expression written in terms of 
the su(2) generators J°, J 1 * 1 acting in the auxiliary space (we will also use the notation J l k ). 
It is clear from u(l) invariance that any such O has a definite grading under J°, namely: 

[J°, O] = s(0) O . (C.l) 

If s(0) 7^ 0, then its trace will vanish for any complex value of the spin. So let us consider 
the case s(0) = 0. In this case O can be rewritten as 

0^F(J°,J 2 ), (C.2) 

where F is some function. 

C.l Twist regular izat ion and </> — > limit 

Let us consider 

Tr^t (e- i<t)j0 o) - Tr^, i (e'^^o) , (C.3) 
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where the trace is taken over the infinite dimensional Verma module. We will show that 
this quantity is finite in the — > limit. 

To proceed it is useful to realize the highest weight representation explicitly in a Fock 
space. The latter is an infinite module 7T Z for any I. It thus contains the information about 
the spin in our realizations of the generators, see (I3.3ip and H3.32J) . As just discussed, we 
will take O as in ( IC.2D . Then (IC.3D becomes 



oo 



J2 ( e -^ n -^F(n-j,j(j + 1)) - e ~^ n+ ^F(n + j + 1, j(j + 1))) . (C.4) 



n=0 



Introducing f{n — u) = F(n—j,j(j + l)), where v = j + h, and we suppress the dependence 
on the Casimir j(j + 1), (1C.4J) can be rewritten as 



oo 



e _i 5 J2 {e~ iHn ' u] f(n - v) - e~^ n+v) f{n + v)) . (C.5) 

n=0 

Consider the case f(s) = s k for some integer k. Then (1C.5I) reads 

-" K) E (— - — ') ^ «^ K)* (^Sf) ■ <-) 

This expression is finite in the — )■ limit for any fe, z/. This concludes the proof of the 
finiteness of Tj(z; 0) 

C.2 Trace functional 

We now want to show that the procedure of evaluating traces over complex spin of the 
previous section is the same as the one used in [5U] and termed "trace functional" . It has 
the property thato 

rrv functional t Jo _ smn (t V) (H7\ 

ir ^ 6 = ^rnh|~- [L - 7} 

This formula plays the role of a generating function for the traces of (Jo) fe - To show the 



equivalence with the definition of the previous section it is enough to consider (1C.3J) with 

O = 1 and — i(p = t. This gives 

smh{tv) 



sinhf ' {Ci 



and completes the proof of equivalence. 



27 For convenience of the reader the notation used here is related to the one used in [30] in the following 
way: x = 2v, H = 2 Jo, z = |. See [3D] for the complete set of defining properties of (|C.7p . 
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C.3 Other approaches 

We will now show that there is a third way to define traces over complex spin, used in [27]. 
It is equivalent to the previous two procedures. Now the definition of the trace is 

Tr ^° 0(3[) = j dfi( % 7) Otfjp) (7 r( 7 ^ - , (C.9) 



where 

^(7,7) =e~ T7 ] [ , k , 7-7 = 7 1 7 1 + 7 2 72 . (CIO) 

■*- ■*- /7T7 



2m 
fc=i 



As discussed before only the zero weight part of 0(J l k ) will contribute to the trace, cf. f lC.lj) . 
To make contact with the previous analysis we compute 



t 



Tr Pronko e tJ 0= I dll^^etV W^* 



,1 9 .,2 a 



(7-7) 



-\2i/-l 



r(2i/) 



l f ( t i 2 - N 2 " _1 

^(7,7) e2 7 7l + e 2 7 72 



r(2i/ 

sinh(t z/ 
sinh I 



(c.11; 



This last equation shows the equivalence with the other two definitions of trace. 

D su(2) covariance properties 

In this appendix we will study some su(2) transformation properties of Q±. Thanks to the 
identity 

e is- ( 1 a+ \ = g _ 7a + / 1 a+ \ g7a + , . 

la^_+a + a" \ a~ z_ + a+ a~ / v ' ; 



where 



^ = ( ' \ ) , (D.2) 



the trace of a monodromy matrix built from 



1 a+ 

a - 2_ + a + a" 



(D.3) 



will satisfy the equation 

e 7S ^ Q_(*_) = Q-0-) . (D.4) 

This is a formal equation because for the trace to give a finite result one needs to introduce 
a regulator e~ l( ^ a a , and this operator does not commute with e ±7a . Still one can write 

e 75 - Q_(^_; <P) = Q-(z_; 0) (I + 0{(f>)) . (D.5) 

33 



In a similar way one can show that 

Q+(^ + ; 0) e^ 5 - = Q + (^+; 0) (I + 0((f>)) . (D.6) 

On the other hand Q± are u(l) invariant for any value of the twist, namely 

e 75 - Q±(*±; 0) e^ s ?<* = Q ± (z ± ; 0) . (D.7) 

Using the previous identities one can show that 

e**- Q±(4; 0) Qi x (4; 0) e~^ tot = Q±(4; 0) Q± x (4; 0) (i + o±(0) • (d.s) 

The proof is trivial using that Q±(£ x , 0) and Q^ 1 ^ 2 , 0) commute. So the su(2) invariance 
in the — )• limit follows. 

E More on sl(2) transformation properties of Wron- 
skian 

Let us investigate how Q zo (z; 0) and P zo { z ', 0) transform under rotations of Q±. One has 



One then finds 



Q 1-+W6 (q :»' With «*-^=l. (E.l) 



P^; 0) -)• P Zo (z; 0) , Q 2o (z; 0) -)■ Q 20 (z; 0) + P zo (z; 0) ©(? s \z ; 0) . (E.2) 



The explicit form of O is not important here, just notice that O vanishes for 

a fi\ (a \ A* /A _/0 ia~ l \ /OiWa 0\ 

7 5 J "" l a^ 1 i ' U(jJ _ Ua J^'UoJlOar 1 ]' 



(E.3) 



which is the reason why equations (14.121) . (I4.13P cannot be inverted for Q±. This freedom 
corresponds to the normalization of Q± and to the discrete operation of exchanging Q + 
with Q_. 

F Numerics 

In this appendix we want to present more details on the method used in section 14.51 for 
finding the two-cut root distributions for both the Q and the P polynomials. We will 
perform it in two steps. Firstly, we will find all roots of the polynomial Q(u) using the 
method presented in [31]. Then, using the information about the polynomial Q(u) we will 
find P(u) from the generalized Wronskian relation (14.141) . We will end this appendix with 
a few figures presenting different aspects of the problem. 
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According to [3T], in order to find roots of the polynomial Q(u) for given M and L ~ M, 
we have to take as a starting point a slightly modified system. Let us fix the number of 
excitations M, and choose the length of the spin chain L' ~^> L such that 4j -C |. One 
finds that the approximate root distribution in this case is given in terms of the roots Zk 
of Hermite polynomials as 

u k = — ^- (l + ^v 7 ^ + 0(1)) . (F.l) 

Here, Z\~ are the solutions of H{zk) = 0, where by H we denoted a Hermite polynomial. To 
simplify our considerations we will investigate only the simplest two-cut solution with mode 
numbers given by n& = 1, for k = 1, . . . , =£, and n^ = — 1, for fc = 4f + 1, . . . , M. Knowing 
the approximate solution, we can use it as a starting point for the Newton method. This 
way we can find the solution of (12. 7p with desired precision for any M and much larger 
L'. Unfortunately, this method does not work for the configurations close to half-filling 
because then the approximation (IF.ip is not good enough, and the Newton algorithm no 
longer converges. However, to overcome this difficulty, we can treat the root distribution as 
a function of the length, while keeping the magnon number M fixed. Assuming in addition 
that the root distributions for slightly differing L's should be numerically close, we can 
use the solution found for a given length as initial data for Bethe equations with smaller 
length. This way we get a sequence of configurations with fixed M and decreasing L, 
corresponding to increasing filling ^. This procedure is sufficient to produce the solutions 
of Bethe equations with excitation numbers very close to half-filling. In particular, we can 
get the configuration presented in figured] in the main text. 

In order to get the root configuration of the P-polynomials we will use (j4.14p with j = 
and in the </>—>■ limit. At the left hand side we then have only u L , while at the right hand 

side we can substitute Q(u) = Y\i 2 =i{u 2 — uf) . Here we denoted by U{ all roots with positive 
real parts found at the previous stage. We assume the polynomial P(u) to be of the form 
P{u) — (u — Uq) J2i=o a i u% wri ere we factorized out [u — Uq) in order to use the freedom 
of P{u) which we discussed in the main text. We can indeed always construct a solution of 
the untwisted Wronskian in the form P'(u) = P{u) + a Q{u) by choosing a such that any 
given Mo become a root of P' . This way we can find a one-parameter family of solutions 
numerated by Uq. Now for given u we substitute the polynomial P(u) of the mentioned 
form into the Wronskian relation. This way we will end up with the set of L — M + 1 
linear equations, where the coefficients aj, i — 0, . . . , L — M will play the role of unknowns. 
The solution of this linear equations is easy to find. It leaves us only with the problem of 
finding the roots of P(u), which we can perform numerically. Using the methods presented 
above we are able to find configurations for chosen values of L and M below half-filling. 
Additionally, we may observe the behaviour of the solution for different Uq, see also [53] . 
In the main text we presented only the Uq = case. In figure [2] we illustrate the root 
configurations also for other values. 
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Figure 2: Root distribution of Q{u) (purple squares) and P(u) (blue dots) with various u 
for L = 60 and M = 16. (a) -S — 0, (b) u = 1, (c) u = 4. Reading this figure from the 
bottom one can see that when |w | is decreased the two condensates shown in (c) start to 
approach each other (b) and eventually assemble on one line (a). It is interesting to note 
that for even length L the two condensates pass through each other without "touching" 
(like in the figure), while for odd L they "scatter", thereby producing double roots on the 
imaginary axis for u = 0. 



36 



References 

[1] L. N. Lipatov, "Evolution equations in QCD", in "Perspectives in Hadronic Physics," 
Proceedings of the Conference, ICTP, Trieste, Italy, 12-16 May 1997, eds. S. Boffi, 
C. Ciofi Degli Atti and M. Giannini, World Scientific (Singapore, 1998). 

[2] V. M. Braun, S. E. Derkachov and A. N. Manashov, "Integrability of three-particle 
evolution equations in QCD," Phys. Rev. Lett. 81 (1998) 2020, hep-ph/9805225. • 
V. M. Braun, S. E. Derkachov, G. P. Korchemsky and A. N. Manashov, "Baryon 
distribution amplitudes in QCD," Nucl. Phys. B 553 (1999) 355, hep-ph/9902375. 
• A. V. Belitsky, "Fine structure of spectrum of twist-three operators in QCD," 
Phys. Lett. B 453 (1999) 59, hep-ph/9902361. • A. V. Belitsky, "Renormalization 
of twist-three operators and integrable lattice models, " Nucl. Phys. B 574 (2000) 407, 
hep-ph/9907420. 

[3] J. A. Minahan and K. Zarembo, "The B ethe- Ansatz for Af = 4 super Yang-Mills," 
JHEP 0303, 013 (2003), hep-th/0212208. • N. Beisert, C. Kristjansen and M. Stau- 
dacher, "The dilatation operator of M = 4 Super Yang-Mills theory," Nucl. Phys. B 
664 (2003) 131, hep-th/0303060. • N. Beisert and M. Staudacher, The M = 4 SYM 
Integrable Super Spin Chain," Nucl. Phys. B 670, 439 (2003), hep-th/0307042. 

[4] I. Bena, J. Polchinski and R. Roiban, Hidden symmetries of the AdSs x S 5 superstring, " 
Phys. Rev. D 69 (2004) 046002, hep-th/0305116. 

[5] N. Beisert and M. Staudacher, "Long-range PSU(2,2\A) Bethe ansaetze for gauge 
theory and strings," Nucl. Phys. B 727 (2005) 1, hep-th/0504190. • N. Beisert, 
B. Eden and M. Staudacher, " Trans cendentality and crossing," J. Stat. Mech. 0701 
(2007) P021, hep-th/0610251. 

[6] G. Arutyunov and S. Frolov, "String hypothesis for the AdS§ x S 5 mirror," JHEP 
0903 (2009) 152, arXiv : 0901 . 1417 [hep-th] . • N. Gromov, V. Kazakov and 
P. Vieira, "Integrability for the Full Spectrum of Planar AdS/CFT, " arXiv : 0901 . 3753 
[hep-th] . • D. Bombardelli, D. Fioravanti and R. Tateo, "Thermodynamic Bethe 
Ansatz for planar AdS/CFT: a proposal," arXiv: 0902. 3930 [hep-th]. • N. Gro- 
mov, V. Kazakov, A. Kozak and P. Vieira, "Integrability for the Full Spectrum of 
Planar AdS/CFT II," arXiv: 0902. 4458 [hep-th]. • G. Arutyunov and S. Frolov, 
"Thermodynamic Bethe Ansatz for the AdS 5 x S 5 Mirror Model," JHEP 0905, 068 
(2009) , arXiv : 0903 .0141 [hep-th] . 

[7] R. J. Baxter, "Partition function of the eight-vertex lattice model," Annals Phys. 70 
(1972) 193 [Annals Phys. 281 (2000) 187]. 

[8] V. V. Bazhanov, S. L. Lukyanov and A. B. Zamolodchikov, "Integrable Structure of 
Conformal Field Theory II. Q-operator and DDV equation," Commun. Math. Phys. 
190 (1997) 247, hep-th/9604044. 

37 



[9] V. V. Bazhanov, S. L. Lukyanov and A. B. Zamolodchikov, "Integrable structure of 
conformal field theory. Ill: The Yang-Baxter relation, " Commun. Math. Phys. 200 
(1999) 297, hep-th/9805008. 

[10] V. V. Bazhanov, A. N. Hibberd and S. M. Khoroshkin, "Integrable structure of W(3) 
conformal field theory, quantum Boussinesq theory and boundary affine Toda theory, " 
Nucl. Phys. B 622 (2002) 475, hep-th/0105177. 

[11] V. V. Bazhanov and Z. Tsuboi, "Baxter's Q-operators for supersymmetric spin 
chains," Nucl. Phys. B 805 (2008) 451, 0805.4274 [hep-th] . 

[12] V. Kojima, "Baxter's Q- operator for the W-algebra W N ," J. Phys. A. 41 (2008) 
355206(16)pp, 0803.3505 [nlin] . 

[13] H. Bethe, "Zur Theorie der Metalle. I. Eigenwerte und Eigenfunktionen der linearen 
Atomkette", Zeitschrift fur Physik 71 (1931) 205. 

[14] A. N. Kirillov and N. A. Liskova, "Completeness of Bethe 's states for generalized 
XXZ model," J. Phys. A 30 (1997) 1209, hep-th/9403107. • A. N. Kirillov and 
N. A. Liskova, "Completeness of Bethe 's states for generalized XXZ model, II," 
hep-th/9607012. 

[15] R. J. Baxter, "Completeness of the Bethe ansatz for the six and eight-vertex models, " 
J. Statist. Phys. 108 (2002) 1, cond-mat/0111188. 

[16] K. Fabricius and B. M. McCoy, "Bethe 's equation is incomplete for the XXZ model at 
roots of unity," J. Statist. Phys. 103 (2001) 647, cond-mat/0009279. 

[17] V. Tarasov, A. Varchenko, "Completeness of Bethe vectors and difference equations 
with regular singular points',' Int. Math. Res. Notices (1995) 637-669, • E. Mukhin, 
V. Tarasov and A. Varchenko, "Bethe Algebra of Homogeneous XXX Heisenberg Model 
Has Simple Spectrum", 0706.0688 [math]. 

[18] G. P. Pronko and Yu. G. Stroganov, "Bethe Equations "on the Wrong Side of 
Equator" ", J. Phys. A 32 (1999) 2333, hep-th/9808153. 

[19] N. Beisert and R. Roiban, "Beauty and the twist: The Bethe Ansatz for twisted Af — 4 
SYM," JHEP 0508 (2005) 039, hep-th/0505187. 

[20] R. J. Baxter, "One-dimensional anisotropic Heisenberg chain. " Annals Phys. 70 
(1972) 323-337. 

[21] L. D. Faddeev, "How Algebraic Bethe Ansatz works for integrable model, " 
hep-th/9605187. 

[22] V. V. Bazhanov, S. L. Lukyanov and A. B. Zamolodchikov, "Integrable structure 
of conformal field theory, quantum KdV theory and thermodynamic Bethe ansatz," 
Commun. Math. Phys. 177, 381 (1996), 9412229. 

38 



[23] C. Korff, "A Q-operator for the twisted XXX model," J. Phys. A : Math. Gen. 39 
(2006) 3203, math-ph/0511022. 

[24] A. N. Kirillov and N. Y. Reshetikhin, "Exact solution of the integrable XXZ Heisenberg 
model with arbitrary spin. I. The ground state and the excitation spectrum, " J. Phys. 
A 20 (1987) 1565. 

[25] S. E. Derkachov and A. N. Manashov, "Factorization of the transfer matrices for 
the quantum sl(2) spin chains and Baxter equation, " J. Phys. A 39 (2006) 4147, 
nlin/0512047. 

[26] V. B. Kuznetsov, M. Salerno and E. K. Sklyanin, "Quantum Backlund transformation 
for the integrable DST model," J. Phys. A 33 (2000) 171 solv-int/9908002. 

[27] A. E. Kovalsky and G. P. Pronko, "Baxter Q-operators for the integrable discrete self- 
trapping chain, " Theoretical and Mathematical Physics, vol. 142, issue 2, pp. 259-269, 
nlin/0203030. 

[28] In preparation. 

[29] G. P. Pronko, "On the Baxter's Q operator for the XXX spin chain, " Commun. Math. 
Phys. 212 (2000) 687, hep-th/9908179. 

[30] H. Boos, M. Jimbo, T. Miwa, F. Smirnov and Y. Takeyama, "A recursion formula for 
the correlation functions of an inhomogeneous XXX model, ", hep-th/0405044. 

[31] T. Bargheer, N. Beisert and N. Gromov, "Quantum Stability for the Heisenberg Fer- 
romagnet, " New J. Phys. 10 (2008) 103023, 0804 . 0324 [hep-th] . 

[32] N. Beisert, J. A. Minahan, M. Staudacher and K. Zarembo, "Stringing spins and 
spinning strings," JHEP 0309 (2003) 010, hep-th/0306139. 

[33] N. Gromov and P. Vieira, "Complete 1-loop test of AdS/CFT," JHEP 0804 (2008) 
046, 0709 . 3487 [hep-th] . 

[34] R. J. Baxter, "Eight vertex model in lattice statistics and one- dimensional anisotropic 
Heisenberg chain. 2. Equivalence to a generalized ice-type lattice model, " Annals Phys. 
76, 25 (1973). 

[35] R. J. Baxter,, "Exactly Solved Models in Statistical Mechanics,", Academic Press 
(1982). 

[36] K. Fabricius and B. M. McCoy, "New Developments in the Eight Vertex Model," J. 
Statist. Phys. Ill (2003) 323, cond-mat/0207177. • K. Fabricius and B. M. Mc- 
Coy, "New Developments in the Eight Vertex Model II. Chains of odd length, " 
cond-mat/0410113. • K. Fabricius, "A new Q-matrix in the Eight- Vertex Model," 
J. Phys. A 40 (2007) 4075, cond-mat/0610481. 

39 



[37] V. V. Bazhanov and V. V. Mangazeev, "Eight-vertex model and non-stationary Lame 
equation," J. Phys. A 38 (2005) L145, hep-th/0411094. • V. V. Bazhanov and 
V. V. Mangazeev, "Analytic theory of the eight-vertex model, " Nucl. Phys. B 775 
(2007) 225, hep-th/0609153. 

[38] L. N. Lipatov, "High-energy asymptotics of multicolor QCD and two-dimensional 
conformal field theories," Phys. Lett. B 309 (1993) 394. • L. N. Lipatov, "High- 
energy asymptotics of multicolor QCD and exactly solvable lattice models, " arXiv:hep- 
th/9311037. • L. N. Lipatov, "Asymptotic behavior of multicolor QCD at high energies 
in connection with exactly solvable spin models," JETP Lett. 59 (1994) 596, [Pisma 
Zh. Eksp. Teor. Fiz. 59 (1994) 571]. 

[39] L. D. Faddeev and G. P. Korchemsky, "High-energy QCD as a completely integrable 
model," Phys. Lett. B 342 (1995) 311, hep-th/9404173. 

[40] E. K. Sklyanin "Backlund transformations and Baxter's Q-operator,", nlin/0009009. 

[41] M. Gaudin and V. Pasquier, "The periodic Toda chain and a matrix generalization of 
the Bessel function's recursion relations," J. Phys. A 25 (1992) 5243. 

[42] A. G. Choudhury, A. R. Chowdhury, "Canonical and noncanonical variables, Baxter's 
Q-operator and the XXX model," Chaos, Solitons and Fractals 16 (2003) 79. 

[43] V. V. Bazhanov and Yu. G. Stroganov, "Chiral Potts model as a descendant of the six 
vertex model," J. Statist. Phys. 59 (1990) 799. 

[44] G. von Gehlen and V. Rittenberg, "L^ Symmetric Quantum Chains With An Infinite 
Set Of Conserved Charges And Z N Zero Modes," Nucl. Phys. B 257 (1985) 351. 

[45] H. Au-Yang, B. M. McCoy, J. H. H. perk, S. Tang and M. L. Yan, "Commuting 
transfer matrices in the chiral Potts models: Solutions of Star triangle equations with 
genus > 1," Phys. Lett. A 123 (1987) 219. 

[46] R. J. Baxter, J. H. H. Perk and H. Au-Yang, "New Solutions Of The Star Triangle 
Relations For The Chiral Potts Model," Phys. Lett. A 128 (1988) 138. 

[47] S. E. Derkachov, "Baxter's Q-operator for the homogeneous XXX spin chain", J. 
Phys. A 32 (1999) 5299, solv-int/9902015. • S. E. Derkachov, "Factorization 
of the R-matrix. F, J. Math. Sci. 143 (2007) 2773, math/0503396. • S. E. Derkachov, 
"Factorization of the R-matrix. IF, J. Math. Sci. 143 (2007) 2791, math/0503410. • 
S. E. Derkachov, "Factorization of R-matrix and Baxter's Q-operator," J. Math. Sci. 
151 (2008) 2848, math/0507252. 

[48] G. P. Pronko, "On Baxter Q-operators for Toda Chain,", nlin/0003002. 

[49] N. Gromov, V. Kazakov and Z. Tsuboi, "PSU(2, 2|4) Character of Quasiclassical 
AdS/CFT, " 1002 . 3981 [hep-th] . 

40 



[50] D. V. Bykov and S. Frolov, "Giant magnons in TsT-transformed AdS§ x S 5 ," JHEP 
0807 (2008) 071, 0805.1070 [hep-th] . 

[51] A. Voros, "Spectral zeta functions," in " Zeta functions in geometry," eds. Kurokawa, 
N. and Sunada, T. Adv. Stud, in Pure Math. 21, pp. 327-358, Kinokuniya, Tokyo, 
1992. 

[52] P. Dorey and R. Tateo, "Anharmonic oscillators, the thermodynamic Bethe ansatz, 
and nonlinear integral equations," J. Phys. A 32, L419 (1999), hep-th/9812211. 

[53] V. V. Bazhanov, S. L. Lukyanov and A. B. Zamolodchikov, "Spectral determinants 
for Schroedinger equation and Q-operators of conformal field theory, " J. Statist. Phys. 
102, 567 (2001), hep-th/9812247. 

[54] D. Gaiotto, G. W. Moore and A. Neitzke, "Wall- crossing, Hitchin Systems, and the 
WKB Approximation, " 0907 . 3987 [hep-th] . 

[55] L. F. Alday, D. Gaiotto and J. Maldacena, "Thermodynamic Bubble Ansatz," 
0911 .4708 [hep-th] . • L. F. Alday, J. Maldacena, A. Sever and P. Vieira, "Y-system 
for Scattering Amplitudes, " 1002 . 2459 [hep-th] . 

[56] S. L. Lukyanov and A. B. Zamolodchikov, "Quantum Sine (h)- Gordon Model and Clas- 
sical Integrable Equations, " 1003 . 5333 [math-ph] . 

[57] V. V. Bazhanov, S. L. Lukyanov and A. M. Tsvelik, "Analytical results for the Coqblin- 
Schrieffer model with generalized magnetic fields," Phys. Rev. B 68, 094427 (2003), 
cond-mat/0305237. 

[58] N. Andrei, K. Furuya and J. H. Lowenstein, "Solution Of The Kondo Problem, " Rev. 
Mod. Phys. 55 (1983) 331. • A. M. Tsvelik and P. B. Wiegmann, "Solution of the 
Kondo problem for an orbital singlet," J. Phys. C15, 1707 (1982); Adv. in Phys. 32, 
453 (1983). 



41 



